#ifndef __GLOBAL__PARAMETERS__
#define __GLOBAL__PARAMETERS__

#include <cmath>
#include <array>
#include <cstdlib>
#define triangle 3
#define quadrilateral 4

int element_shape = triangle;// triangle or quadrilateral

int degree_of_basis_function = 1; // 1 or 2

double RHS(double x, double y) // right hand side function in Poisson equation
{
    return 2.0 * sin(x) * sin(y);
}

double dirichlet(double x, double y) // boundry function of Dirichlet problem
{
    return sin(x) * sin(y);
} 

#else
//do nothing

#endif
